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We study the Casimir effect for a 3-d system of ideal Bose gas in a slab geometry with Dirichlet 
boundary condition. We calculate the temperature(T) dependence of the Casimir force below and 
above the Bose-Einstein condensation temperature(T c ). At T < T c the Casimir force vanishes as 
[tjt-] 3 '' 2 . For T > T c it weakly depends on temperature. For T ^> T c it vanishes exponentially. At 
finite temperatures this force for thermalized photons in between two plates has a classical expression 
which is independent of h. At finite temperatures the Casimir force for our system depends on h. 
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Vacuum fluctuation of electromagnetic field would 
cause an attractive force between two closely spaced 
parallel conducting plates. This phenomenon is called 
Casimir effect and this force is called Casimir force 
0, EL HI- In the original paper [l| the Casimir force at 
zero temperature (T — 0) was defined as 
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where E(L) is the ground state energy (i.e. the vacuum 
energy) of the electromagnetic field in between the two 
conducting plates separated at a distance L. This force 
has been measured experimentally However, Casimir 
effect can be generalized [5J for any range of temperature 
and for any dielectric substance between two dielectric 
plates. It has also been generalized for thermodynamical 
systems [f|. Casimir force for this kind of systems has 
recently been measured 0- At finite temperature T, the 
definition of Casimir force is generalized as [1, H Ell 
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where £It(L) is the grand potential of the system con- 
fined between two plates separated at a distance L. 

We consider the Casimir effect for a thermodynamical 
system of Bose gas between two infinite slabs. Geometry 
of the system on which some external boundary condi- 
tion can be imposed is responsible for the Casimir effect. 
Thermalized photons (massless bosons) in between two 
conducting plates of area A at temperature T gives rise 
to the Casimir pressure [H El El El O 
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where k is the Boltzmann constant, c is the velocity of 
light and L is the separation of the parallel plates. At 
T —> 0, Casimir pressure becomes — 2 ^ ^ and it is only 
the vacuum fluctuation which contributes to the Casimir 
pressure. At high temperature i.e. for — > 0, the 
Casimir force for photon gas goes as L~ 3 and has a purely 
classical expression independent of h. 



Let us consider a Bose-gas is confined between two in- 
finitely large square shaped hard plates of area A. The 
plates are along x-y plane and they are separated along 
z- axis by a distance L. For the slab geometry, \f~A ^> L. 
We consider that our system is in thermodynamic equi- 
librium with its surroundings at temperature T. At this 
temperature the thermal de Broglie wavelength of a sin- 
gle particle of mass m is A = 



Th 2 



, , , . In the thermo- 
For this system the single par- 



dynamic limit, 7 < 1. 

tide energy is e(p x ,p y ,j) = + f£ + , where 

p x and p y are the momentum along x-axis and y-axix 

respectively and j = 1,2,3, However in the thermo- 
dynamic limit the single particle energy can be written as 

<(pT-.Py:Pz) - — - J ^ 

along z-axis. 



., + 7T~ + where p z is the momentum 
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Considering the thermodynamic limit the total number 
of thermally excited particles can be written as 



N T = 
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where [i is the chemical potential and V is the volume 
of the system. Bose condensation temperature (T c ) is 
defined as a temperature where all the particles are ther- 
mally excited and below that temperature a macroscopic 
number of particles come to the ground state [H, El EH • 
At T < T c the chemical potential goes to the ground state 
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energy. So 
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Let us now introduce the finite size correction. The 
ground state energy of our system is [g = 2m l 2 ] • The 
average no. of particles with energy e Px , Py ,j is given by 

where /j,' = (fj, — g) < for 



bosons. At and below the condensate temperature fi' — » 
0. For this bosonic system we have the grand potential 

fl = fl(A,L,T, //) 
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Replacing j by (j' + 1), we recast the above equation as 
n(u,L,T, fj,') 

Adp x dp y 
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where A = \J 2 ^kT ■ Integrating over p x and p y we get 
il(A,L,T,n') 
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Since ^ -C 1 , higher order terms of the above series would 
not contribute significantly. From Euler-Maclaurin sum- 



mation formula we convert the summation over j' to in- 
tegration. So from equation(9) we have 
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Collecting the leading terms from the above equation(lO) 
we can write 

il(A,L,T,n') 

~ 2ttH 2 ^ i 2 ^Xi 1 / 2 ~ 2 + 
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where z = e M / fcT is the fugacity and gi(z) — z + ~c + 

|j- + .... is the Bose- Einstein function. From the above 
equation we get the total number of particles as 
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In the thermodynamic limit of a system, as T < T c , 
z — > 1. For a finite system this can not happen, other- 
wise the correction terms in the above expression would 
be infinite. Instead at T > T c , z ~ 1. Taking only 
the first correction term in the eqn.(12) we have gi(z) — 
-ln(l-z) = [N'(T)g3(z) - N((3/2)]± = -InAz, where 

N'(T) = 4^[^C(3/2)] and Az = 1 - z is a small 
change in the fugacity at T >T C . Now putting z = 1 in 
the expression of Az, we get Az = e - AAr C(3/2)L/A^ w here 
A7V = N'(T) — N. We see that in the thermodynamic 
limit (L — > oo) Az = and when L is finite such that 
L/X > 1 we have z ~ 1 at T > T c . 

Let us now calculate the Casimir force. At T < T c we 
put /x' — » or z — ► 1. So from eqn.(ll) we have 

^.l.T.O). - ^=[^(5/2) -> 

+ |C(3/2)j] (13) 
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Here the first term of eqn. (13) is 
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It is the bulk term of the grand potential. From our con- 
sideration of thermodynamic limit ^ = constant. So 
Qt(oo) = O;,. The second term of eqn. (13) is (fl s ) = 

A ^a ro [|C(2)]. It is the surface term of the grand po- 
tential. The third term of eqn. (13) is the Casimir term 
of the grand potential. We call it the Casimir potential. 
Now putting N = 4*^2^(3/2) in eqn.(13) we find 
the Casimir potential as 
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Putting A = \J ^mkT m equation(15) we have 



¥ r T i3/2 71-2 h2 
[ tJ 2mL 2 



(15) 



(16) 



1 

0.8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 



FIG. 1: Casimir force(j p c (T, L)) versus temperature(T) plot. 
F C (T, L) is in units of and T is in units of T c . Equa- 

tion (18) corresponds to the range if- < 1. Equation (20) 
corresponds to the range 1 < ^- < 1.1. 



From eqn. (2) and (14) we have 
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F C (T,L) = -—Cl c 
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For T < T c , from eqn. (16) and (17) we have the expres- 
sion of Casimir force as 
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This expression for the Casimir force shows that, at finite 
temperatures the force depends on H. 

Above the condensation temperature \i' < or z < 1. 
However, for T > T c , z ~ 1. So at T > T c , from the 
eqn. (11) and eqn. (5) with trivial manipulation we get 
the Casimir potential as 
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In the above equation for T > T c , we put z = 1. 
However as ^- increases z decreases. So, for T > T c 
the Casimir potential weakly depends on temperature. 
Putting T = T c + AT in equation (19) and from the 
definition of Casimir force we have 
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where < ^ < 1. 



For T > T c , the Casimir force 



weakly depends on temperature. 



Let us now calculate the Casimir force at T ^> T c . 
At these temperatures z«l So we can approximately 
write gi{z) « z + §<-. From the first term of eqn. (12) we 

have gz{z) = N ^TmL ~ z + z /( 2 ^)- For this ran § e 01 
temperatures we can write e^l kT = ze~ 7Tx2 / L ' 2 rts z. For 
convenience, we replace fi' by u and recast equation (7) 
as 
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use the 
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formula £~=-oo e"™" = 
e * From the above equa- 

tion we choose the Casimir potential as f2 c = 

A(kT) 2 m 



where we 
1 v-^oo 



i=l 2^7 = 1 



For 



oo, m 
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the expression of above Casimir potential we can put 
e n/kT — z ^ 1 an d can take j — 1 anc | j = 1 as the 
leading term to contribute in the Casimir potential. So, 
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for T>T C the Casimir potential is 
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equation (22), for T > T c we have the Casimir force as 
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Now we see that in the classical limit (T 3> T c ) the 
Casimir force vanishes as e~ kT . 

The changes of Casimir force with temperature for the 
range < T < T c and for the range T > T c is shown in 
FIG 1. 

That vacuum fluctuation causes Casimir force is well 
known [l], 0]. Critical fluctuation also causes Casimir 



force@,i. The Casimir force calculated here is neither 
due to vacuum fluctuation nor due to critical fluctuation. 
It is due to quantum fluctuation. This fluctuation is asso- 
ciated with the commutator algebra of position and mo- 
mentum operator as well as with the commutator algebra 
of bosonic annihilation operator(eti) and creation opera- 
tor (a|) such that [aj,aj] = 8^j, where i,j represent the 
single particle energy states. At T <C T c almost all the 
particles come down to the ground state. The quantum 
fluctuation dies out due to the macroscopic occupation of 
particles in a single state. That is why the Casimir force 
dies out at T < T c . At T > T c the Bose- Einstein statis- 
tics becomes classical Maxwell-Boltzmann statistics and 
thermal fluctuation dominates over the quantum fluc- 
tuation. For this reason the Casimir force dies out at 
T T c . Below T c the reduction of thermodynamic 
Casimir force with T 3 / 2 law is the signature of Bose- 
Einstein condensation. 
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